Electron capture rates in a plasma 
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A new general expression is derived for nuclear electron capture rates within dense plasmas. Its 
qualitative nature leads us to question some widely accepted assumptions about how to calculate 
the effects of the plasma on the rates. A perturbative evaluation, though not directly applicable to 
the strongly interacting case, appears to bear out these suspicions. 
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INTRODUCTION 

Calculation of electron capture rates in plasmas is of 
importance in many astrophysical systems. The influ- 
ence of the surrounding plasma on the capture rates is 
completely understood for the case of a weakly coupled 
plasma in which the electrons are not degenerate, as long 
as the energy release in the reaction, Q, is sufficiently 
large. These conditions are well fulfilled in the solar 
core, where the capture rate of electrons on ^Be has been 
studied in detail. Here the weak-screening Salpeter fac- 
tor of a multiplicative factor of exp(— c^Zku/T) k, 
(1 - e'^ZKo/T) gives a 16% rate reduction 0, d,!!. ^ 
These calculations of reduction factors are all really just 
calculations of a screening-induced reduction of the ex- 
pectation of the electron density at the position of the 
ion. The methods of the above references are not appli- 
cable to cases with: a) strong plasma coupling F; or b), 
small energy release Q in relation to other energies in the 
problem (or negative Q). Indeed, as we explicitly show 
below, the plasma-induced change of the electron density 
at the position of the ion is no longer the determining fac- 
tor in the rate when condition b) prevails. Also, of the 
approaches in the above references, only that of rcf. [4;] 
can incorporate high electron degeneracy^ 

In the parameter region relevant to explosions/collapse 
of 0-Ne-Mg stellar cores, for example, both conditions a) 
and b) apply and the electrons are very degenerate. In- 
deed Q is here negative for important capture reactions; 
they become enabled when the Fermi energy increases 
nearly to the value —Q. In ref. @, addressing these 
systems, we find a simple statement of a procedure that 
appears to be widely followed 0], 0: 0: [1] ; namely, 
to calculate the capture rate for an individual electron 
as in vacuum, except replacing the vacuum Q value for 
the decay by an effective value Q + Ei — E2 where Ei 
and E2 are the respective energies of interaction of the 



initial and final ions in the plasma. There exist more or 
less standard expressions for these energies derived from 
classical simulations. For the purpose of our present dis- 
cussion we designate this correction, with energies calcu- 
lated as for a classical plasma, as the "ionic energy correc- 
tion." In applications, it appears that this correction has 
been applied sometimes with accompanying initial state 
Coulomb effects involving the electrons, which we loosely 
characterize as "screening" , and sometimes without. 

The present work begins with the derivation of a new 
expression for the rate of electron capture, one that we 
think should be the basis of future calculation of rates 
in the presence either of strong degeneracy, or of strong 
plasma coupling. Just looking at the form of the result, 
we believe, should be enough to make one profoundly 
skeptical about assumptions being made in the literature. 
We then carry out a systematic perturbative evaluation 
in the weakly coupled case, for cases of arbitrary electron 
degeneracy, retaining terms of relative order c'^Zkd/T in 
the rates (or of order "Salpeter" as we shall henceforth 
designate them). The results appear to be in complete 
disagreement with the assumptions that lead to the ionic 
energy correction, had they been applied to a weakly 
coupled problem. 

We use the notation kd for the screening wave number 
throughout. In the presence of some degree of degeneracy 
its square is given by [loj . 



= 47r 



d 



dp 



(1) 



where = eZs is the charge of ionic species s, the 
number density of that species, and f{Ep) the Fermi dis- 
tribution. In the non-degenerate limit this gives the De- 
bye wave number, and in the strongly degenerate limit, 
the electronic part would give the square of the Thomas- 
Fermi screening wave number. 



^ The question of whether or not the exponentiation of the pertur- 
bative result is justified is numerically irrelevant in this case. 

^ although we know of no very interesting set of conditions that 
are highly degenerate without also embodying condition a) or 
condition b). 



GENERAL FORMULATION 

Our general methods follow those introduced in ref. 
[loj |. We take the electron capture interaction Hamilto- 
nian for the process of a non-relativistic electron plus ion. 
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Ja, going into neutrino plus ion to be given by a zero 
range Fermi type coupling (all the effects of Coulomb in- 
teractions in what follows being the same for the Gamow- 
Teller case), 



Hw^ J dr[K{r,t)+K\r,t)l 



(2) 



where 



K{r, t) = ge-'«Vj (r, t)^e(j, t)^l{r, t)^^, t), (3) 

and Q is the energy release in the reaction. There are 
suppressed spin indices contracted with each other on the 
electron and neutrino fields. Here the V'{a,b} cire nonrela- 
tivistic quantum fields that describe creation or annihila- 
tion of the respective ions. All fields are in a Heisenberg 
picture with respect to the complete Hamiltonian. The 
ionic fields could be Fermi or Bose, but we shall remain in 
a domain of temperature and density in which deviation 
from Boltzmann statistics of the ions is inconsequential. 

We are interested in the time rate of change of the 
electron density induced by ([2]). Since the medium is 
translationally invariant we can choose to evaluate this 
time derivative at at point r = and choose time to be 
zero as well. Directly from the Heisenberg equations we 
obtain the rate of change, w, of the electron density at 
the origin, ne(0, 0) = ^1(0, 0)V'e(0, 0), 

w=(^he{0,0))^^-i J dr{J[ne(O,0),if(r,0)]^^ 

= z([i^(O,0)-ift(O,0)])^. (4) 

The notation (...)^ indicates the thermal average in 
the medium, such that for an operator, O, we have 
(0)/3 = ZpiTr[Oexp(-/3[ff + Hw - ^J-eNe])] where Zp 
the partition function, and where H + Hw is the com- 
plete Hamiltonian. We wish to calculate the rate to low- 
est non-vanishing order; i.e. to second order in the weak 
coupling parameter g. So it is it clear that one power of 
g must come from the weak interaction Hw within the 
statistical factor. Thus we now must consider the lin- 
ear response of the medium average of the operator 
[K{Q) - K^Q)] to this perturbation, 

^ = ~^ f_^dt i\dv)l^[n,{Q,Q),Hw{v]) ^ = - jl^dt 
X / dr([K{0, 0) - K^O, 0)] , [i^ (r, t) + K'<{r, t)] 



(5) 

where now the thermal average in the medium is to be 
calculated using H alone. H conserves electron number, 
whence, 

([i^(O,0),i^(r,t)]) =([i^t(o,o),ift(i.,t)]^ ^0. (6) 



Using in addition the space-time translational invariance 
of the medium and the antisymmetry of the commutator. 



7. t(o, 0), i^(r, t)])^ = -([if (0, 0), K^i-v, -t)])^, (7) 



we can write the rate as 



w 



dt I dY([K{0,Q),K\Y,t) 



/3 



(8) 



When we take the medium to contain no ions of type 6, 
so that there is no reverse reaction, we can omit the first 
term in the commutator in ([8]), which, after inserting ([3|). 
gives. 



w^g^ I dt dre~^Q'Ul{r,t)i^i{r,t) 

J —OO J 

X V',(r, t)Mr, t)i^l{0, 0)Vj (0, 0)^,(0, 0)Ve(0, 0)^ 

d^p^ 
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dt / dreg's* 



(27r)6 



X (V4(r, t)^l{r, t)M^, i)v4 (0, 0)Va(O, O)^e(0, 0) 



(9) 



For almost any purposes the thermal motions and re- 
coil of the reacting ions are ignorable. We distinguish the 
ion la on which an electron is captured and the capture 
product lb from the other ions in the plasma, taking a 
single /„ fixed at the origin, replacing r = in all of 
the explicit fields in (|9]) and adding an all-over factor of 
density of ions of type a, n^. The products of the field 
operators ipa,b that occur in ([9]) now serve only to change 
one ion into the other. We substitute 



^\{Q,t)MO,t)^e^"'\a){b\e 



iHt 



(10) 



where the time dependence refiects the fact that the 
plasma interactions of the two ions are difi^erent. We 
obtain, 



27r2 



dte-'^' / dp^p^e'^"'M{t), (11) 



where 



Tr 



M{t) = Zp^ X 
(ar^«e*^VI(O,0)|a)(6|e-'^*|6)(a|Ve(O,0)|a) 



(12) 



The ionic states |a),|5), in ([T2|) serve to determine the 
coupling to the fixed ion in the interaction Hamiltonian. 
The trace in (fT2|) is over all of the other coordinates. We 
write 



M{t) 



Zp^Tr 



e-(/3-»t)^fa^t(o,0)e-''^''Ve(O,0) 



(13) 
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where 



Zp = Tr e-^ff^+^/^e^. 



(14) 



Here Ha is the Hamihonian for the initial system (with 
Hw turned off, since our expressions are aheady of order 
g^) and Hf, that for the final system, not counting the 
constant Q. The space of states in which the trace in 
([T3|) is to be evaluated is just what we have defined as 
the "plasma" ; all electrons, and all ions except the ions a 
and b. To restate the definitions of Ha and iJf,, they are 
the Hamiltonians for the complete plasma in the presence 
of the respective fixed Coulomb fields of the nuclei, a, b. 

The result (fT3|) for the function M{t), in conjunction 
with the rate expression ()11|) is our fundamental result. 
We show the details of the perturbation evaluation in the 
appendices, retaining terms of order Salpeter (e^ZKu). 
This calculation appears to be very inefficient, given the 
simplicity of the result, but we have found no other path 
to the answer. We state the results in terms of two sub- 
sidiary functions, Mo{t) and M^(t), 



Mo{t) = 



and 



dpi dp 

(27r)3 



MAit) 



iE„t 



d^p 
(2^ 



fiEp)e 



iE„t 



(15) 



d^p 



[/(i?p)]'e'^''* 



(16) 



(27r)3 ' 

where f(p) is the Fermi distribution, 

/(ii;p) = [l + e^(^--^=)]-i. (17) 
The kernel of the rate expression (|lip is now given by, 

Mit) = Mo{t) -eeanoP^Moit) - MaII)), (18) 
and the rate by, 

w = 23/277^2^3/2 / dE £;i/2 (Q + E) 

Jmin(0,-Q) 

x'f{E)-eeanDf3[f{E)~f{E)]\. (19) 

We note that to the order that we are working, the re- 
sult p9)) agrees exactly with eq. 1.21 of ref. [4| except 
for the additional term, E in the factor Q -|- £^ in our 
result, and the inclusion of all vacuum Coulomb effects 
in ref. (In the present work we just stated that 

vacuum Coulomb effects could be treated additively.) In 
the vacuum case the extra E term just puts in corrected 
kinematics. But, as we see more explicitly in the next 
section, once we make this trivial modification, the rate 
is no longer determined by the electron density at the 
location of the ion. Refs. 0] and Q calculate only the 
perturbation of the density. Thus it is of some interest 



that this modification is screened by the same factor (un- 
der the E integral) as the Q term, but it is perhaps not 
unexpected. What is more of note is not what our result 
contains, but what it does not contain. We discuss this 
at some length in the next section. 



DISCUSSION 

We return to the basic result for the kernel of the cap- 
ture rate formula, 



M{t) = Zp^Tr 



■/3+»t)Jfo^t(o,0)e-'^''*Ve(0,0) 



(20) 



If we were to replace both factors exp{iHat) and 
exTp{—iHf)t) by unity, i. e. evaluate at t — 0, then M 
would be precisely the electron density at the position of 
the ion. Then, in accord with the remarks at the end of 
the last section, we get the standard screening results, 
where Q + E in is replaced by Q. In our mechanics 
this replacement comes from the fact that Hamiltonians, 
Ha and Hb both contain electron kinetic energy terms, 
and that there is an electron creation operator standing 
between the factors containing them; thence producing 
an exp[—iEt]. But Ha and H}, also contain pieces in 
which the initial and final ions interact with the plasma. 
So in parallel fashion, why do we not see another addition 
to Q in the results, namely AQ = AEa ~ AEh, due to the 
ionic interaction shifts in the plasma? After all, one can 
infer energy shifts for the ions at the Salpeter level. For 
an ion with charge we have AEa — ^a'^D /2 so that 
the difference between the energy shifts of the initial and 
final ions in our case would be eCaUD- But there is no 
such term in the result (fT9|) . 

This is our first criticism of the "ionic energy shift" 
ansatz. It should have manifested itself at our level of 
approximation and it did not. A reader might comment, 
"I look at p9)) . say in the non-degenerate limit, and I see 
just Salpeter screening, an effect that pertains only to 
physics in the initial state. But we know that the ener- 
gies of these ions are shifted differently by the plasma in- 
teractions, and when these shifts are of the order Q, they 
must matter a lot. So the author cannot have put in all 
of the physics." But if one wades through the calculations 
in the appendix, one finds that such terms come and go; 
they all cancel in the end. What is left out of the ansatz, 
[initial state effects -I- final state effects], is the fact that 
many of the individual terms in the development, even 
to this order, involve one Coulomb interaction with the 
initial ion and one with the final one. These terms can- 
cel almost all of what comes from the ansatz. We see no 
valid reason to believe that the ionic energy shift correc- 
tion by itself, as it is applied in refs. [l], [3|: 01; will 
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serve us any better in the strongly coupled case than in 
our case in which we can really calculate. 

In addition to this ionic energy shift, there exists a sec- 
ond plausible correction recommended in some papers 
dealing with extremely dense and degenerate systems 
, [1] , namely an electron energy shift resulting from the 
screening of the potential seen by the electron. The qual- 
itative description of how this screening should affect the 
rates is now nothing like the description in weak screen- 
ing theory; in the latter case the effects are entirely due 
to the change in electron density at the ionic position 
and are not directly connected to the energy spectrum of 
the electrons. But when we address a degenerate case in 
which Q < and we get captures only as the Fermi en- 
ergy approaches —Q, it is the energy distribution rather 
than the electron density that would seem to be the crit- 
ical feature. 

To discuss this suggested effect, we first note that in 
the derivation of our perturbative results, given in the ap- 
pendix, we did need to calculate a change in the electron 
chemical potential, (5/ie, Our primary calculation was for 
fixed chemical potential and we thereafter made a cor- 
rection to get the results for fixed electron density. The 
outcome is that we obtain the Salpeter correction with 
no explicit Sfi appearance, with the stipulation that the 
chemical potential to be used is obtained from the num- 
ber density by doing the integral rig — J dp{2n)^^ f {Ep) 
and then solving for fi^, as a function of rie- It appears 
to us that this gives the statement of results in the best 
form for applications, and there is no explicit chemical 
potential shift in the result. 

Thus we find reasons to question both of the above 
corrections; they do not contribute in the weakly coupled 
case, and we see no reason for them to be the dominant 
correction in the strongly coupled case. That is not to 
say that the physical reasoning behind them is incorrect; 
other corrections are equally important. 

Turning to specific applications we first focus our at- 
tention on ref. which is a comprehensive treatment 
of electron capture in the infall phase of a type II su- 
pernova event, with a full range of nuclear species, tem- 
peratures and densities. Here the screening corrections 
are a relatively minor feature of the work; still they give 
effects ranging from a 25% decrease in rates for a density 
of lQ^^gc~^, temperature .75MeV (with plasma param- 
eter r « 1 ) to a factor of two at a density of 10^^, 
T = 1.5MeV, r w 3 . Our suggestion here is that these 
corrections be taken as an estimate of the magnitude of 
the unknown error in either direction due to plasma ef- 
fects. If this turns out to create important uncertainties 
in the outcome of the collapse, then it will be time to 
worry. 

In another venue, the accretion of matter onto the sur- 
face of a neutron star, the problems could be much more 
serious. Here, at a depth where nuclear reactions stimu- 
lating transient explosive phenomena can occur, it may 



be that electron capture on a particular nuclear species 
is the stimulator and the rates will be extremely density 



and temperature dependent [13|. In this environ- 

ment, the plasma couplng F will be much stronger than 
those cited in the previous example, and we can expect 
the effects on rates to be larger, and harder to estimate, 
even in magnitude. Furthermore, it is the temperature 
dependence more than the absolute rate that rules in 
most explosion calculations. Since the plasma correc- 
tions, whatever else, will be strongly temperature depen- 
dent, we could say that, knowing nothing about them, 
we equally know nothing about the effective exponent in 
the temperature dependence of the rate. 

There is a body of literature that presents unorthodox 
results on the effects of the plasma on reaction rates, even 
in the case of the weakly coupled solar plasma. These re- 
sults, if correct, would have significant impact on solar 
models. We wish it to be clear that the work of the 
present paper does not support any of these claims. One 
can find a critique in ref. lj| as well as a list of some 



of the articles in this category, though there are a more 
recent examples that go down the same path [l^ . There 
are, in addition, a few papers of the same genre address- 
ing electron capture [16|, 17|, lisl . The last of these 
engendered a response in ref. [l9|, one with which we 
completely agree. But, that said, the reply in [l^ was 
more in the vein of "We did a correct calculation and 
had no such terms; therefore your calculation is wrong." , 
rather than tracking down specific mistakes. The same 
could be said of the arguments given in ref. Q . Thus it 
may be worth adding my own conclusion about one thing 
that has gone wrong in at least some contributions of this 
genre. Let me characterize as an "S-matrix approach" 
one in which one starts the correction of a two body 
reaction rate by putting in some more bodies from the 
plasma in the calculation, working out the multi-body re- 
action rate in some perturbative way and then doing the 
thermal average (typically) over the momentum states of 
the plasma particles. Looking at the actual mechanics in 
these papers, it appears to me that the assumptions are 
tantamount to entirely leaving out Coulomb forces in the 
factor exp[—/3Ha] in (IT3|) . It is a mistake to do this. We 
could put the conclusion qualitatively as: A treatment 
strictly from a multi-particle S-matrix, followed by aver- 
aging over thermal distributions, misses essential physics 
coming from the fact that the initial particles sit in an 
interacting medium 



THINKING ABOUT THE STRONGLY 
COUPLED AND VERY DEGENERATE CASE. 

Ionic energy shift 

Returning to the basic equation (pil)) that determines 
the kernel of the capture rate formula, we go back and 
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think qualitatively about the ionic energy shift. We first 
state the case for the shift assumption, and then the 
counter-argument. We take only electron kinetic energies 
and ionic potential energies in Ha.b, the electron potential 
energies not being relevent in the extremely degenerate 
system, at least not to the effect that we are consider- 
ing. Now we look first at the Hi, where there is a term 
efc^/(0), the energy of ion b in the fields of all of the other 
ions. If we forget about where b came from, and when 
it came into being, and consider a problem in which it 
has been at r = all along, we can in principle calculate 
the potential of the plasma ions at the position of b; it is 
of course a function of the charge = eZt- In classical 
simulations the energy shift is given by AEh = T g{Th) 
where 17(2;) is of order x and = Z^^'^e^3^^/'^(47rne)^/'^. 
We can calculate the corresponding energy shift for ini- 
tial ion a. From the time dependent exponentials in (pil)) 
we would infer that these corrections change the effective 
Q value by an amount AEa — AEb- 

But the actual term in the interaction, say. Ha, that 
induces an energy shift for ion a is ea0(O). First, through 
this term, the plasma redistributes itself around the ion, 
creating a potential, which then shifts the energy of the 
ion through another action of Ha- In a perturbative cal- 
culation the plasma is thereupon returned to its original 
state. If we let Ha act only once, it will transfer a bit of 
momentum from ion a to the plasma. Of course the whole 
expression is a trace, so everything has to be returned to 
its initial value. But that bit of momentum can be re- 
turned to the plasma particles after the electron capture, 
through the interaction H},. In the weakly coupled case 
there is substantial cancellation between these categories 
of effects. We cannot argue that in the strongly coupled 
case they would cancel in any particular way. But we 
would not trust a result that takes one part and discards 
the other 

Ab initio numerical calculations? 

We take a box containing Nj ions, not including a 
or 5, and N^. electrons. The initial and final ions a, b 
that participate in the capture process are at the point 
r = 0, and manifest themselves only in the potentials Ha 
and Hi,. In one respect we are at an advantage here, 
in comparison with a parallel treatment of the nuclear 
fusion problem, because we probably do not care about 
regions in the space of the ions, Ri, that are classically 
inaccessible, or nearly unaccessible at the temperature of 
the medium. Thus we can plausibly take the ions to be 
absolutely classical, and entirely omit the ionic kinetic 
energy terms from Ha and H/,. We take a basis set of 
Ne electrons at positions ri, r2....rAr^ , and Nj ions at 
positions Ri, R2....R7Vj- Using the cyclic invariance of 
the trace in ([^0]) to move -0(0,0) from the last position 
to the first, and putting in the position basis we have. 



M{t) = Zp^J dRi, ..(IRnj, ..dr2, --drN^ J dr'2, --dr' x 
({0,r2, .rArJ, Ri, .Rw,|e-('3~»*)«''|{0, r'2, .r'^JRi, -R^,) 
x({r'2, ..r'jvj, Ri, ..RAr,|e-**^''|{r2, ..r^J, Ri, ..Rat,). 



We have here shifted from the grand canonical ensemble 
to the canonical ensemble, with a definite number of elec- 
trons. The { } around the electron basis states signify 
a total anti-symmetrization. The ionic part requires no 
symmetrization. Note that the intermediate integration 
is over electronic positions only. We have taken the clas- 
sical limit for the ions, in which the ionic kinetic energy 
term is omitted and the ionic coordinates within Ha and 
Hb remain at their external values. 

Could we in principle use he trick of ref. [iot . 
exp{-l3H) = [exp{-l3H/Ns)]^^ taking Ns is sufficiently 
large so that we can make a lowest order perturbative ex- 
pansion in the Coulomb potentials in each individual fac- 
tor? Each factor would then be at a sufficiently high tem- 
perature, NsT, to validate this lowest order expansion. 
This would be at best an expensive calculation because 
A^i, Ne and Nj all have to be fairly large, and between 
each of the A^i factors in the product one must integrate 
over the full manifold of some r[...r'j^ , not to mention 
the enormous anti-symmetrizations required. And the 
final kiss of death would be: because of the real time ex- 
ponentials, we would have rapidly oscillating integrands. 

But perhaps instead, we could start by taking a very 
degenerate plasma. The electron sea is so stiff that it 
does nothing, except that one gets captured, creating an 
electron hole in the final state. Should we be worried 
about the effects of the (ionically) irregular plasma on 
what that hole does? Perhaps somewhat symbolically, 
we can write a much simpler equation, 

M{t) = Zp^ f dRi, ..dRNj (Ri, ..Rat, |e-(/5-**)^» |Ri, .Rn,) 
X (r 0, Ri, ..R^v, |e-**^' |r = 0, Ri, ..Rat,). 

(22) 

Is it possible to calculate the matrix elements for the hole 
to go from the origin back to the origin in time t under 
the action of Hi, , for a big set of irregular distributions of 
ions ? If one had to know the answer; the eikonal method 
might serve. 

This work was supported in part by NSF grant PHY- 
0455918. 



APPENDIX. PERTURBATION CALCULATIONS 

To demonstrate the methods we divide the respec- 
tive Hamiltonians, as Ha^t = Hq + H^j^ , in which 
the interaction term contains all Coulomb interactions 
of the electrons in the plasma and the interactions of 
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the distinguished ions la and /& with both electrons and 
plasma ions. Coulomb interactions among the ions in the 
plasma are retained in Hq. We denote the charges of 
the electron and two ions respectively as e, Bq, Cf,, where 
66 = Ba — 1 obtaining, 



where 



Hi 



and 



Hi 



fdr 



Jdr 



e[0/(r) + i0e(r)]ne(r) + e/n/(r)(/)(r) 

+ea</'(0), (23) 



e[0/(r) + i0e(r)]ne(r) 

+660(0). 



e/n/(r)(/)(r) 



(24) 



Here (/)e, 4>i are the respective electric potentials pro- 
duced by the electrons and </> = + </>£ is the total 
potential of the plasma. We have not included the ba- 
sic Coulomb interactions between electrons and la and 
lb- But if we repeated the calculations to come using 
Coulomb wave- functions instead of plane waves we would 
not find any medium-dependent terms that are of leading 
order, e^Ku. There is of course a conventional Coulomb 
correction in the absence of the plasma; perturbatively 
it begins at order e^, and for solar core capture in ^Be is 
larger than the order c^kd screening corrections. At the 
level of the present paper this vacuum contribution can 
just be added to the answer; Coulomb wave-functions are 
not needed for calculating the rest of the corrections. 

We introduce the interaction picture through the iden- 
tities. 



and, 



where, 



and 



^-H^ifi~^t) ^ e-^°(^-'*)l)+(-i/3, t) , (25) 



e-^""* = n_{0,t)e-'""\ (26) 



l]i+)(-*/3,t)=exp * / dt'Hiit') , (27) 



H^^(t) = e'Ho{t+^0)HI^-^Hoit+^fS)_ (28) 

In ((27)) the subscript, stands for time ordering of the 
integrals. Similarly we have. 



nl \0,t) =exp -i dt'Hl{t') 
L Jo 



(29) 



^ In our terminology, the distinguished ions la and are not part 
of the "plasma" . Since Fermi statistics prevents us from distin- 
guishing an electron; we include all electrons in the "plasma" . 



(30) 



and now the integral is anti-time-ordered. Substituting 
in (fT3l) we obtain 



M{t) = Zp^Ti e-f^»('5-**)VKO,0)6-^-f^^*Ve(0,0) 



x^e(O,0) 



e-^«'3f]i+V*/3,0V4(O,0) 



(31) 



We expand the n factors in powers of the coupling, 
retaining only terms of order e^, (or e^, eca, 6^, etc). 
Schematically, each will have powers of e coming from 
the thermal expectation of a product of two ion-electric- 
potential (/>(r) operators; these are determined by the 
Hamiltonian Hq. If this potential-potential correlator is 
itself expanded in powers of e, the expansion begins with 
terms of order , so that the rate corrections would be of 
order but the terms have an infrared singularity from 
the long range Coulomb force, and when that is properly 
regulated the entire correction to the rate will be of order 
e^. In the present paper we pursue only these terms, 
the neglected terms will go as power or higher. For 
this purpose we shall need the part of the field correl- 
lation function that is the most singular as goes to 
zero. 



(32) 



which is time independent. 
We now expand fii"^^ 



dhHiih) f ' dt2Hi{h), 

il3 J-i/3 



(33) 



and Q 



(-) 

6 ' 



f](-)(o,t) = i + i / dhiiiih) 

t nti 

dtiHiiti) / dt2Hl{t2). 
^0 



(34) 

Before getting bogged down in the evaluation of a mul- 
tiplicity of terms we think for a moment about the gen- 
eral form of the expressions generated by the expansion. 
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Leaving out most of the arguments of the functions and 
not concerning ourselves with orders of factors, we note 
that the form is the expectation value, 



time independent, there is in effect no time dependence. 
Putting dSS]) into dSS]) and dM]), and then putting the 
results into (|3T|l and using ([32| we obtain. 



1 



(35) 



Now we state some results of peering into the term by 
term outcome: 

1. ) The electron operators within ?/'e,V'l each get 
paired with an electron operator in an rig or term ear- 
lier in the expression. The terms in which they eat each 
other are "disconnected" , i.e. their contribution gets can- 
celed by a perturbation in the partition function in the 
denominator of (1311) . in conjunction with an unperturbed 
numerator. 

2. ) Terms where both i/j and ^/i^ pair with the operators 
from the same big square bracket in (j35|) are not infrared 
divergent in the absence of screening and will therefore 
not contribute in leading order. 

3. ) Bearing in mind 2.) above, when a single ip or i/j^ 
is to be paired with an operator from the product. 



i y"dr(/)e(r)ne(r) = J dr J d.r'\ 



^ne(r)ne(r'). 



(36) 



then the pairing will be, in effect, the same when it is with 
an operator in rig or when it is with one in (j)e. Thus if 
we add a stipulation to the calculation of ((35|) that the 
explicit and ip'^ are never to be paired with (j>e in this 
kind of term, we can compensate by multiplying the 0e 's 
in (1351) by 2. To put it more simply, we can replace Ha^b 
of ^ and dMl) by. 



Hi= J dr[-e0(r)ne(r)]+e,0(O), 



H^, = I drhe0(r)ne(r)] +660(0), 



(37) 



with the above stipulation on computation. These oper- 
ators H^fj{t) in the interaction picture are now given in 
terms of the electron creation and annihilation operators 
a^p),a{p) by. 



ea,b 



(27r) 



3 "p+k"P 



i2ny 



(38) 



where t ^ t + i/3 for the incoming system a and t ^ t for 
the final system. In the the small k limit we can set t — 
in the final bracket in ([38|) . We do keep an infinitesimal 
k in the indices of the operators Cp+k, as we shall see 
below. Since the infrared part of the correlator is also 



M{t) = Mo{t) -4tt 



(27r)3 k^{k^ + Kjy) 







2 )Mi{t) -(t + ^/3)^Af2(^) + jM^it) 



(39) 



where 



Mi{t) = Limfc^o 



dp4 a 



P4-k'^P4 



Ga - e y dps a^''' 
dpidpe'^"*"^ 



P3+k"P3 



a^^flp), (40) 



M2(t) = Limfc^o 



^P3+k'^P3 



dpial^ 



66 



dp4 al 



P4-k°'P4 



dpe^^-*ap), (41) 



A'hit) = Limfe^o{^ / c'Piaj,! ^6 - e / dpa Op 
X 66 - e/dp4a|,^_k«P4 \dp. 
and A/q is the unperturbed part, 



P3+k'^P.3 



■P/' 



(42) 



Mo(i) = J dpidp(al^ape'' 



Ept 



(43) 



Here f{p) the Fermi distribution, f{Ep) — (1 + 
g0(-Ep-/ie)-)-i g^jj^ ^j^g electron chemical potential is de- 
termined from the electron density to zeroth order in e 
by Ue = / dp{2n)^^ f {p) . When we calculate perturba- 
tive corrections we must explicitly include a correction 
that readjusts the relation between chemical potential 
and density, which we shall do below. 

The thermal averages in (pUf - P^ are readily calcu- 
lated. A number of terms are disconnected and canceled 
by perturbative corrections to the partition function: a) 
all terms of order in which p = pi, i.e. in which op- 
erators in V'(O) a-i'e paired with those in V'^(O); b) some 
more that are discussed below. The rest of the terms can 
be expressed in terms of Mq of the function M^, 



2 lEpt 



(44) 



and a time independent constant, R, 



{ea - e fMoit) + (2eae - e^)MA{t) + e^RMo{t) , (45) 
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where i? is a time independent constant that is propor- 
tional to the system volume and that will be canceled by 
a corresponding perturbation of the partition function. 

We reconstitute the complete time dependent function 
M{t) which when inserted in (|lip gives the rate. 



M{t) 



'■Kofi 



Moit) 



+ -{2eea - e')KDpMA{t) - -Kope" RM^it)]. (46) 

Now we calculate the changes of the same order to the 
partition function, where the system is governed by Ha^ 

I-ifjdh f_]pdt2 X [ea(j}iO,ti) + ef dr(j){r,ti)ne{r,ti)] 
x[ea</)(0,i2)+e/ drV(r',i2)ne(r',t2)]])^. 

(47) 

The calculation follows the lines outlined above, but 
it is simpler. We obtain a contribution of order that 
is the same as the order term in ()46p . but without 
the factor of Mo{t). Therefore in conjunction with the 
1 X Mo(t) term in ([31]) . we see that the partition function 
perturbation removes the e^KD/3Mo(t) term from (|46|. 
The other non- vanishing term coming from (j47p is of or- 
der and it exactly removes the final term, proportional 
to R, from (gH). 

Next we must recognize that in our calculation the 
chemical potentials were held fixed while we calculated 
perturbations to the rates. But the relation between 
number density and chemical potential has corrections of 
exactly the same order as the rates. In the next section 
of the appendix we calculate the changes in this relation, 
and we find that the re-expression of the results in terms 
of the corrected densities removes the remaining terms 
of order from (|46l) ^. leaving the final result for the 
corrections of leading order as simply, 

M{t) = Mo{t) - e ea^Dli(Mo{t) - MA{t)) . (48) 

Keeping rie fixed 

In sec. 4 we used a relation between the change in the 
electron number density, Srif,, when we turn on Coulomb 
interactions in the medium keeping the chemical poten- 
tial fixed, and the chemical potential change that 



* This mirrors the situation in the case of a perturbative calcu- 
lation of plasma corrections ot fusion rates, as discussed in ref. 



is then required to restore the original electron density. 
Then this inserted as a correction into the Fermi 
distribution in the unperturbed rate function ^o(^) gen- 
erates a change that must be added to the terms shown 
in ([i^ . in order that the complete result give the change 
in rate induced by the Coulomb forces while keeping the 
number density, rather than the chemical potential, con- 
stant. We have 



(49) 



where the operator f2^(0, — z/3) is expanded in exactly the 
form shown in psp . except that we omit the interaction 
of the ion with the plasma in the calculation; that is we 
take Ca — Eb — Only the second order in the expansion 
contributes, and we see that the answer is given in terms 
functions that we have already exhibited, now evaluated 
at i = 0, 

ne = 4"^ + ^e2^D/?[Afo(0) - M^(0)], (50) 

This gives a density shift that is compensated by a chem- 
ical potential shift. 



A 1 2 



(51) 



We must compensate with a shift 5^e in the function 
Mo(t). Expanding, this gives a shift in the function that 
determines the rate. 



5eM{t) 



1 



(52) 
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